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Abstract 

The defect induced chemical expansion in acceptor-doped barium zirconate is investigated using density-functional theory (DFT) 
calculations. The two defect species involved in the hydration reaction, the -b2 charged oxygen vacancy and the proton interstitial 
forming a hydroxide ion, are considered both as free defects and in association with the dopants Y, In, Sc and Ga. The defect induced 
strain tensor X is introduced, which provides a natural generalisation of the ordinary chemical expansion to three dimensions and to 
anisotropic distortions. Both the addition of a vacancy and a proton cause anisotropic distortions and a net contraction of the lattice, 
indicating that both the vacancy and the hydroxide ion are smaller than the oxygen ion. The contraction is considerably larger for 
the vacancy and the net effect in hydration, when a vacancy is filled and two protons are added, is an expansion, consistent with 
the experimental findings. The effect of the dopants on the chemical expansion in hydration is found to be quite small, even if it is 
assumed that both the vacancy and the proton are fully associated with a dopant atom in the lattice. 
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l. Introduction 

Over the past three decades acceptor-doped perovskites have 
been studied extensively for the use as proton conducting elec¬ 
trolytes in intermediate temperature fuel cells, electrolyzers, 
etc. HI El. One of the most promising materials is BaZrOs, 
which combines high bulk proton conductivity together with 
chemical stability towards CO and CO 2 a. 

Protons are incorporated into the perovskite structure through 
the hydration of oxygen vacancies, which are formed due to 
the acceptor doping in order to maintain charge neutrality. In 
Kroger-Vink notation the hydration reads as 

H20(g) + v^-+0^^20H^ (1) 

It has been seen experimentally that the hydration of acceptor- 
doped BaZrOs is associated with a volume expansion GliHTl 
that puts the material under mechanical stress and can lead to 
micro-cracking of the dense electrolyte and ultimately deteri¬ 
oration of the cell performance. Chemical expansion due to 
oxidation has been reported for other well-known ionic con¬ 
ductors, such as Ce 02 -j 11113 and Lai_;,Sr;,Coi_yFe 3 , 03_5 d 

m. In these materials the expansion is related to two phenom¬ 
ena: (1) the formation of oxygen vacancies and (2) the reduc¬ 
tion of cations, which change the ionic radius of those ions, 
e.g. Ce^^ ^ Ce^^. The latter effect is not present in BaZrOs 
as the barium and zirconium cations are isovalent. Thus, the 
hydration induced volume expansion should be more directly 
related to the size of oxygen vacancies and proton interstitials 
(hydroxide ions). 
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The size of the hydroxide ion is well known in terms of the 
Shannon ionic radius, and it is slightly smaller than the oxygen 
ion E2 However, the size of the charged oxygen vacancy is 
not as well defined. Suggestions have been made that the for¬ 
mation of a vacancy should lead to an expansion due to the 
electrostatic repulsion between the cations located closest to 
the vacancy nmia. On the other hand, studies of Ce02 and 
Zr 02 indicate that the vacancy is smaller than the oxygen ion 

Eiiiiiiiia. 

In this paper we study the size of the -\-2 charged oxygen 
vacancy and the proton interstitial in acceptor-doped BaZrOs 
using density-functional theory (DFT). The size of the two de¬ 
fects is determined in terms of the defect induced strain tensor. 
This quantity does not only predict the lattice expansion, "the 
size ” of the defect, but also anisotropic expansions, ’’the shape ” 
of the defect. Furthermore, the implications of the results on the 
chemical expansion due to hydration are discussed both in the 
dilute limit of freely moving defects and as trapped defects in 
association with different dopants. 

2. Theoretical formalism 

2.7. Chemical expansion coefficient 

To quantify the chemical expansion induced in a crystal 
upon formation of point defects we start from a classical ther¬ 
modynamics point of view cni. We consider the volume, V, to 
depend on temperature, T, pressure, P, and defect concentra¬ 
tion Xk of type k 

V=V(T,PAxk}) 

The unit for the defect concentration, Xk, is the number of de¬ 
fects per primitive unit cell (one formula unit BaZrOs) and the 
volume of one primitive unit cell is denoted In analogy 
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with the thermal expansion coefficient a = V~^{dVldT)p^Xk 
compressibility k = -V~^{dVldP)T,Xk we define the chemical 
expansion coefficient as 


V [dxkjpj^ 


( 2 ) 


The coefficient Pk corresponds to the chemical volume expan¬ 
sion of the material and in general we have Pk = PkiJ, P, 

The differential for the volume can be written as 





( 3 ) 


To describe the expansion in a solid material including aniso¬ 
tropy, we generalise the concept of a scalar volume and intro¬ 
duce the elastic strain tensor e. The volume expansion, or vol¬ 
ume dilatation, of a material under assumption of small strain 
is given by the trace of the strain tensor e ifT^fTTl 


the displacement gradient can be written as 


(ll;' -1) 


ij 


dui 

drj 


and describes the deformation of the crystal CD. In addition to 
strain, the displacement gradient can also include rotations. Ro¬ 
tations should be eliminated, since they do not alter the internal 
relationships between the atoms and are of no thermodynamic 
significance. In general this can be done by means of a po¬ 
lar decomposition dl, but for small rotations it is sufficient to 
add the transpose and take the mean. Hence the strain tensor is 
given by 


^ij = 


1 

2 [drj dri) 


(5) 


or in matrix notation 


e=i(LLoU(Lo-ifL^)-l (6) 




= Tr(e) 


and the chemical expansion coefficient can thus be written as 


A = Tr 


— ) 


(4) 


If no rotations are present the strain tensor e equals the displace¬ 
ment gradient 

e=(L-Lo)Loi (7) 

This expression is the natural generalisation of the one dimen¬ 
sional strain s = Isala^ to higher dimensions and is similar to 
the one proposed by Centoni et al. da. 


2.2. Strain Tensor 

In order to derive an expression for the strain tensor we 
consider the deformation of a homogeneous anisotropic elastic 
medium. Let 

be the basis vectors defining the undeformed crystal and Lq the 
one-to-one matrix transformation that maps the cartesian coor¬ 
dinates 0 / onto the crystal a^. After the deformation the crystal 
is defined by the new basis vectors 

a', = Le, = LLg^a, 

The location of a point, P = {p\,p 2 ,p?), vn \hQ undeformed 
crystal can then be written as 

i 

and after the deformation 



The displacement of the point is thus 

u = r' - r = ^ p, (a; - a,) = (LLq' - i) r 

i 

where the I is the identity matrix and the quantity (LL"^ - 1^ 
is recognised as the displacement gradient in the Lagrangian 
description of continuum mechanics ITSl . In tensor notation 


2.3. Defect induced strain tensor 

Consider now the effect of a single defect introduced in a 
volume Vq. The defect induced strain tensor A is defined as 

A=-e (8) 

Xd 

where Xd = TldVo is the defect concentration with Tic the vol¬ 
ume of the primitive unit cell. Hence the chemical expansion 
coefficient in Eq. 0 can be written as 

/3 = Tr(A) (9) 

The tensor A provides a natural generalisation of the ordinary 
chemical expansion coefficient to anisotropic deformations. It 
is identical to the T-tensor introduced by Nowick and Berry 
1201 . and directly related to the defect strain tensor A = Tic A 
used by Freedman et al. ED. 

3. Computational method 

Electronic structure calculations were performed using density- 
functional theory (DFT) as implemented in the VASP software 
(version 5.3.3) 122^251 . which uses a plane-wave basis set with 
periodic boundary conditions in all directions. The generalised 
gradient approximation (GGA) PBE was used for the exchange- 
correlation functional and the projector augmented wave method 
(PAW) EH EH to describe the ion-electron interaction. Two su¬ 
percell sizes were used, consisting of 4 x 4 x 4 and 2x2x2 
primitive unit cells, which implies defect concentrations of 1/64 
and 1/8 respectively. The corresponding k-point sampling was 
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2x2x2 and 4x4x4 using the Monkhorst-Pack scheme. The 
energy cutoff was set to 520 eV. To compensate for the elec¬ 
trical charge of the defects a homogeneous background charge 
was added. 

Calculating the strain requires high accuracy. This is due 
to that the strain, which is a small quantity, is obtained as the 
difference between two similar values (see Eq. Q). The elec¬ 
tronic convergence criterion was set to 10“^ eV and the cell was 
allowed to relax in both size and shape until the forces were less 
than 10-3 eV/A. 

The calculations were performed according to the following 
procedure. An ideal bulk BaZrOs supercell was constructed and 
relaxed with resulting supercell parameters Lq. A single defect 
was then introduced into the supercell followed by atomic struc¬ 
ture relaxation where the supercell size and shape were allowed 
to change, giving the defective supercell parameters defined by 
L. The strain tensor e was then obtained from L and Lq using 
Eq. 0. From here the defect induced strain tensor A and the 
chemical expansion coefficient (3 were determined using Eqs. 
0 and 0- For the doped systems we first relaxed the structure 
with the dopant alone, which gave the corresponding chemical 
expansion. A defect was then added followed by further relax¬ 
ation. In this case the doped system served as Lq. 


4. Results 

In this work we consider the two defects associated with 
the hydration reaction in Eq. 0: the +2 charged oxygen va¬ 
cancy (Vq ) and the proton (OHq) ||29l|30||. The computational 
approach makes it possible to study these defects without hav¬ 
ing to consider the dopant atom explicitly by instead adding 
an appropriate homogeneous background charge ED. How¬ 
ever, we have also studied the effect of vacancy-dopant and 
proton-dopant association by explicitly considering four differ¬ 
ent dopants: Y, In, Sc, and Ga. 

The ground state for the defect free system is found to be 
cubic (Pm3m) with lattice constant ao = 4.236 A. This is some¬ 
what larger than the experimental value of 4.19 A but is in 

line with previous GGA results in the literature |[33H36ll . 


4.1. Oxygen vacancy 

First we consider the oxygen vacancy. Fig. illustrates 
the Zr02 and BaO lattice planes, [010] and [100] respectively, 
which intersect at the defect. The [001] Zr02-plane is not ex¬ 
plicitly shown. The ionic displacements in these planes due 
to the formation of an oxygen vacancy are shown in Fig. 
Fig, [^corresponds to distortions in the [010] and [001] planes 


(Zr02) which are equivalent due to symmetry, while Fig. 2b 
corresponds the [100] plane (BaO). 

The ionic displacements are listed in Tab. and are de¬ 
termined with respect to the defect lattice site, the lattice site 
where the defect is created. In general, the displacements can 
be understood from electrostatic considerations. The vacancy, 
having an effective charge of -r2, repels the metal cations and 
attracts the negatively charged oxygen ions. For the nearest 
neighbouring zirconium and barium ions the displacements are 



(a) (b) 


Figure 1: Illustrations of (a) the Zr02-plane [010] and (b) the 
BaO-plane [100]. The red, blue, green and white atoms cor¬ 
respond to oxygen, zirconium, barium and hydrogen, respec¬ 
tively. The added proton shows the defect lattice site as well as 
the orientation of the hydroxide ion. 


0.17 A and 0.10 A respectively, and directed away from the va¬ 
cancy, while the nearest neighbouring oxygen ions are displaced 
0.22 A towards the vacancy. Ions further away from the vacancy 
are displaced as well, although to a lesser extent. 

The displacement of the outermost oxygen ions in Fig. 
indicates the displacement of the boundary of the simulation 
cell and are related to the defect induced strain tensor A. The 
oxygen ions are all displaced towards the vacancy leading to a 
defect induced strain tensor for the oxygen vacancy of 
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-0.038 0 0 ^ 

0 -0.096 0 

^ 0 0 -0.096, 


The components along the diagonal correspond to the usual lin¬ 
ear expansion s = b^aja^ in the different directions divided 
by the defect concentration. All components are negative in¬ 
dicating a contraction of the lattice in all directions. The first 
component is different from the second and third components, 
which are equal due to symmetry, resulting in an anisotropic 
strain. The contraction is smallest in the v-direction, the same 
direction in which the nearest neighbouring zirconium ions are 
being displaced. Since all components of Ay are negative it fol¬ 
lows that the chemical expansion coefficient is negative as well. 
The value we obtain is 


j3y = -0.230 

which indicates that there is a decrease in volume upon for¬ 
mation of a vacancy, or in other words, the oxygen vacancy is 
smaller than an oxygen ion. 

4.2. Proton 

The distortions induced by the proton interstitial are shown 
in Fig. The proton reduces the symmetry with respect to the 
vacancy and consequently the distortions in the [010] and [001] 
planes (Zr 02 ) are no longer equivalent. 

Tab. summarises the proton induced local distortions on 
individual ions. Also in this case the displacements are deter¬ 
mined with respect to the defect lattice site. Similar to the oxy- 
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(b) 

Figure 2: Visualisation of the lattice distortions for an oxygen 
vacancy in (a) the [010] ([001]) plane (Zr 02 ) and (b) the [100] 
plane (BaO). Colours follow the same convention as in Fig. 
and the arrows indicate the displacements larger than 0.05 A in 
Tab.[^ 


Table 1: Displacements of ions relative to the defect lattice site. 
The first three columns state the number of the coordination 
shell (CS) in the ideal cell, the distance to the defect lattice 
site (J) and the atomic species, respectively. The fourth and 
sixth column indicate the coordination number (CN) of each ion 
and the fifth and seventh the length of the displacement vector 
(AJ). The sign indicate whether the radial component of the 
displacement vector is towards (-) or away from (+) the defect. 


CS 

d (ao) 

Atom 

Vacancy 


Proton 

CN 

Ad (A) 

CN 

Ad (A) 

1 

1/2 

Zr 

2 

0.173 (-r) 

2 

0.155 (-r) 

2 

V2/2 

Ba 

4 

0.096 (-r) 

2 

0.198 (-r) 






2 

0.059 (-) 



0 

8 

0.220 (-) 

4 

0.063 (-r) 






2 

0.085 (-) 






2 

0.214 (-) 

3 

1 

0 

2 

0.000 (-r) 

2 

0.159 (-r) 




4 

0.072 (-h) 

1 

0.175 (-r) 






1 

0.107 (-) 






2 

0.093 (-) 

4 

V5/2 

Zr 

8 

0.018 (-) 

4 

0.059 (-r) 






2 

0.040 (-h) 






2 

0.060 (-) 

5 

V6/2 

Ba 

8 

0.023 (-) 

4 

0.065 (-r) 






4 

0.068 (-) 



0 

8 

0.084 (-r) 

4 

0.092 (-r) 






4 

0.053 (-r) 




8 

0.023 (-) 

4 

0.040 (-r) 






4 

0.073 (-) 

6 

V2 

0 

4 

0.026 (-) 

2 

0.027 (-r) 




8 

0.094 (-) 

2 

0.047 (-) 






4 

0.036 (-r) 






2 

0.037 (-) 






2 

0.065 (-) 


defect induced strain tensor for the proton 


gen vacancy the nearest neighbouring zirconium ions are re¬ 
pelled by the proton causing a displacement of 0.16 A. The four 
closest barium ions are also repelled, although two of them are 
substantially more displaced than the others (0.20 A and 0.06 A 
respectively) due to the asymmetric nature of the hydroxide ion 
along the y-direction. For the nearest neighbouring oxygen ions 
the picture is somewhat different compared to the oxygen va¬ 
cancy. The ions in the [001] plane (Fig. are displaced asym¬ 
metrically by 0.21 A and 0.09 A towards the proton, while the 
ions in the [010] plane (Fig. are displaced 0.06 A away. 

The oxygen ions on the cell boundary are also displaced 
differently. There is now an expansion of the supercell in the 
x-direction while in the y and z-directions there is still a con¬ 
traction but not as large as for the vacancy. This leads to the 




^0.009 

0 


0 0 ^ 

-0.066 0 

0 - 0 . 010 , 


where the different components indicate an anisotropic strain in 
all three directions. Since the boundary ions are displaced out¬ 
wards in the x-direction the corresponding strain tensor com¬ 
ponent is positive. Despite this expansion, the contractions in 
the other two directions lead to a negative chemical expansion 
coefficient of 

= -0.066 

Similar to the case for the oxygen vacancy the lattice contracts, 
although the magnitude of the contraction is considerably smaller 
for the proton. This result agrees qualitatively with the Shan¬ 
non radius of the hydroxide ion which is slightly smaller than 
the radius of oxygen ion ifT^ . 
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(C) 

Figure 3: Visualisation of the lattice distortions for a proton 
interstitial in (a) the [001] plane (Zr 02 ), (b) the [010] plane 
(Zr02) and (c) the [100] plane (BaO). Colours follow the same 
convention as in Fig. [J and the arrows indicate the displace¬ 
ments larger than 0.10 A. 


4.3. Acceptor dopants and dopant association 

We have also considered the effect of dopants and dopant 
association. In Tab. [^we show the resulting chemical expan¬ 
sion coefficient from a single dopant yS/)ina2x2x2 supercell 
for four different dopant species. This corresponds to a defect 
concentration of 12.5%, which is within the range of typical 
concentrations (10-20%) for real materials. The cubic structure 
is maintained in all cases. The trend for the chemical expansion 
follows the same trend as the ionic radii by Shannon ifT^ . i.e. 
Y > In > Sc > Ga. 

It is known that there is an attractive interaction both be¬ 
tween dopants and vacancies and between dopants and protons 
Ea and dopant association is expected to be present in the real 
material Wl- The expansion caused by the vacancy and the 
proton may therefore change if these defects are associated with 
a dopant atom. We have determined the chemical expansion for 
a vacancy and a proton with a neighbouring dopant, denoted 
ySv|D and y^HiD respectively, and the result is given in Tab. The 
resulting A tensor is still diagonal for the dopant-vacancy sys¬ 
tem while it becomes slightly distorted with off-diagonal com¬ 
ponents for the dopant-proton systems. 

It is seen from the data in Tab. [2] that the contraction due to 
the formation of the vacancy increases with increasing dopant 
radius. The same trend is observed for the proton interstitial. 
For the largest dopant Y, the distance between the B cation and 
the neighbouring oxygen ion increases in doping from 2.12 A 
for Zr-0 to 2.22 A for Y-0, consistent with experimental data 
|[38l . This difference is less than what is obtained from the 
Shannon radii (0.18 A) and the doped system is under compres¬ 
sive strain close to the dopant. When the vacancy or proton is 
added the strain can partly be reduced and the contraction be¬ 
comes larger compared with the undoped system. 

Table 2: Chemical expansion coefficients for doped and un¬ 
doped BaZrOs in 2 x 2 x 2 supercells. The results for the un¬ 
doped 4x4x4 supercell are included for comparison. The 
second column shows the ionic radius as tabulated by Shannon 
02 and the third column the dopant chemical expansion. The 
ionic radii for the dopants should be compared with the radius 
of 0.72 A of the Zr ion. Column 4 shows the expansion coeffi¬ 
cient for the trapped vacancy relative to the doped supercell and 
column 5 the corresponding coefficient for the trapped proton. 



rs(A) 

ySo 

/?v|D 

^H|D 

undoped^ 



-0.230 

-0.066 

undoped^ 



-0.246 

-0.076 

Y 

0.90 

0.241 

-0.278 

-0.103 

In 

0.80 

0.168 

-0.269 

-0.089 

Sc 

0.745 

0.096 

-0.252 

-0.082 

Ga 

0.62 

-0.000 

-0.242 

-0.074 


^4x4x4 

>’2x2x2 
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5. Discussion 

5.1. Size of the oxygen vacancy 

We find that the lattice contracts when the vacancy is cre¬ 
ated. To quantify this effect the concept of a radius of a spher¬ 
ical vacancy has been introduced for fluorite structured oxides 
(mu da as well as for perovskites Eiiia. In the present 
case we find an anisotropic distorsion which is more correctly 
described by an ellipsoid rather than a sphere. For the change 
of the lattice constant we get in the v-direction ISax = aoXyAxx, 
where Vy is the vacancy concentration and Axx is the first of 
the three diagonal components of the tensor A. In the y and 
z-directions we get the same change Aay = Aa^, with Aay = 

aQXyAyy . 

How to model the distortion of the BaZrOs lattice in a con¬ 
sistent way using the concept of radii is not at all obvious. 
BaZrOs is cubic with a Goldschmidt ratio close to unity. Fol¬ 
lowing Andersson et al. we introduce an effective anion ra¬ 
dius as 

ro,a = I [(3 - Xv)ro + = ro + \xy (rv,a - tq) (10) 

where vq is the radius of the oxygen ion. The subscript a 
denotes the three different directions x, y and z and indicates 
that the effective radius can differ in the different directions. 
In the present case we can write the lattice constant as ax = 
2(rB + fo) in the x-direction and ay = + ^o) foi* the 

y and z-directions, where ta and are the radii of the A- 
site (barium) and B-site (zirconium), respectively. The change 
of the lattice constants becomes Aax = 2xy(ry^x - ^o)/3 and 
Aay = ^^2xy(ry^y “ tq ) / 3 and hence 

rv,y = ro + -^aoAyy ( 12 ) 

Using the Shannon radius ro = 1.40 A da we get for the semi¬ 
principal axes of the ellipsoid: ry ^ = 1.18 A and = 

0.58 A. This can be compared with the value = 1.18 A used 
by Andersson et al. (T) in order to obtain an adequate fit to their 
experimental data. 

The size of the oxygen vacancy has also been the subject 
of interest in studies of other perovskite oxides. Freedman et 
al. (ZD have performed atomistic modelling of SrTiOs using a 
shell-potential model by Akhtar et al. 1^ . They found that the 
vacancy expands in the x-direction while contracts in the other 
two resulting in approximately no change in volume, i.e.jSy ^ 0. 
We have performed calculations for BaZrOs with the same type 
of potential using parameters taken from the work by Stokes 
and Islam ll40l . Like Freedman et al. we found that /3y ^ 0. 

The size of the charged oxygen vacancy has been found to 
be smaller than the oxygen ion also in fluorite structured ox¬ 
ides. Hong and Virkar found the radius of the vacancy to be 
1.164 A and 0.993 A in Ce02 and Zr02 respectively ifT^ . Mar- 
rocchelli et al. found the radius to be 1.169 A and 0.988 A for 
the same materials using both theoretical modelling and exper¬ 
imental techniques Ea. 


5.2. Chemical expansion due to hydration 

With chemical expansion coefficients for the oxygen va¬ 
cancy and the proton interstitial we can now consider the chem¬ 
ical expansion due to hydration. As can be seen from Eq. Q, 
hydration of the material corresponds to replacing one oxygen 
vacancy with two proton interstitials. The chemical expansion 
per incorporated proton is thus given by 


^hydr — 2 ^^ 


(13) 


In Tab. [^we show the result for both supercell sizes. The 
chemical expansion is positive, which corresponds to an expan¬ 
sion of the lattice. The chemical expansion for the two different 
supercell sizes differs by about 4%, indicating an interaction 
effect due to the finite defect concentration. 

The formula in Eq. ( p3] ) neglects the explicit effect from 
the dopants. If association between the dopant and the vacancy 
and/or proton is present the chemical expansion will change. 


However, we found in section 4.3 that the effect is quite small. 
If we assume that the vacancy is associated with a single dopant, 
as a nearest neighbouring defect, the expression in Eq. 03 
should be modified to 


<d.=^H -^AiD (14) 

If also the proton is trapped at a dopant site Eq. 03 should be 
modified to 

Pbydc =/^H|D - jAiD (15) 

In Tab. |3] the results are summarised. The effect from the 
dopants is not pronounced. When only vacancy-dopant associ¬ 
ation is included the chemical expansion increases with increas¬ 
ing dopant radius. If also proton-dopant association is included 
the expansion is rather similar to the undoped system, except 
for the yttrium doped system where the expansion is reduced. 

As mentioned previously in the introduction, lattice expan¬ 
sion due to hydration in BaZrOs has been seen experimentally 
for several different acceptor dopants and dopant concentration 


Table 3: Chemical expansion due to hydration for undoped and 
doped BaZrOs in 2 x 2 x 2 supercells. The results for the un¬ 
doped 4x4x4 supercell are included for comparison. The sec¬ 
ond column shows the expansion according to Eq. The 

third column shows the expansion for trapped vacancies (Eq. 
( p^ ) and the fourth column for trapped vacancies and protons 
(Eq. ((13). 



^hydr 

^hydr 

^hydr 

undoped^ 

undoped^ 

Y 

0.049 

0.047 

0.063 

0.037 

In 


0.059 

0.046 

Sc 


0.050 

0.044 

Ga 


0.045 

0.047 


‘‘4x4x4 

'’2x2x2 
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levels. Based on the measured results we can extract the rel¬ 
ative lattice expansion {/^a/ao) which can be used to calculate 
the chemical expansion coefficient according to 

^hydr — 3 (16) 

■^H 

where xu is the proton concentration. 

Kreuer ^ measured the relative lattice expansion in 15% 
Y-doped BaZrOs and found it to be 0.0040. An almost fully 
hydrated sample with a proton concentration of 0.14 (data ex¬ 
tracted from Fig. 2a) has been reported previously by Kreuer 
et al. BTl . Kinyanjui et al. |[4l studied the hydration of 50% 
In-substituted BaZrOs. The hydration gave rise to an expan¬ 
sion of 0.0107 (data extracted from Fig. 2), although the proton 
concentration was only 0.40 (data extracted from Fig. 3). Fur¬ 
thermore, Hiraiwa et al. O performed a study on 20% Y-doped 
BaZrOs. They found a relative expansion of 0.0019 (data ex¬ 
tracted from Fig. 3) which is a bit lower than the previously dis¬ 
cussed values. They do not present any results related to the de¬ 
gree of hydration although it is quite reasonable to assume that 
the material was not fully hydrated based on the small expan¬ 
sion. This work was followed by a study of 20% doped BaZrOs 
where four additional dopants where considered, namely Sc, 
Sm, Eu and Dy DU. Again, the relative lattice expansions were 
quite small, ranging from 0.0007 to 0.0016 (data extracted from 
Tab. 2). Besides the possibility of an incomplete hydration the 
small expansion values could also be related to the fact the these 
dopants have a tendency to occupy the A-site instead of the 
B-site OHS, which reduces the concentration of oxygen va¬ 
cancies and consequently the maximum proton concentration. 
Additionally, the authors found results indicating that the Dy 
ions could change oxidation state. This is associated with a 
change in the ionic radius, which would affect the chemical 
expansion. Finally, in a recent study by Andersson et al. la 
chemical expansion was observed for several different Y-doped 
samples of BaZrOs with dopant concentrations of 5%, 10% and 
20%. The relative expansion was 0.0015 for both the 5% and 
10% doped samples while it was 0. 0034 for the sample with 
20% yttrium. TGA data presented in the paper indicate that the 
5% doped sample was fully hydrated while the samples with 
10% and 20% dopants were not. 

The chemical expansion coefficient O^hydr) has been calcu¬ 
lated for the discussed experimental results and are summarised 
in Tab. together with the relative lattice expansion and the 
degree of hydration. All these values for y^hydr indicate an ex¬ 
pansion that is consistent with our result in Tab. 

It seems that we may underestimate the expansion based 
on the comparison with the data by Kreuer et al. ^ |4T1 and 
Kinyanjui et al. IH. One uncertainty is the approximation for 
the exchange-correlation functional. The data in this work were 
obtained using the GGA/PBE approximation. However, some 
test calculations using the standard local density approximation 
(LDA) indicate that the results are rather insensitive to the par¬ 
ticular choice of exchange-correlation approximation. Finite 
temperatures and thermal fluctuations may also influence the 
chemical expansion. These effects are not included in this work 


Table 4: Experimental values of the relative lattice expan¬ 
sion coefficient due to hydration (Aalao) and the correspond¬ 
ing chemical expansion coefficient O^hydr) for acceptor-doped 
BaZrOs with various dopant species and concentrations, xa de¬ 
notes the nominal dopant concentration and xu the estimated 
actual proton concentration. 


Dopant 

Xa 

Aajao 

Xh 

^hydr 

Ref. 

Y 

0.15 

0.0040 

0.14 

0.086 

Da ED 

In 

0.50 

0.0107 

0.40 

0.080 

Dl 

Y 

0.20 

0.0019 

<0.20 

> 0.028 

DD 

Sc 

0.20 

0.0014 

<0.20 

> 0.021 

m 

Sm 

0.20 

0.0010 

<0.20 

> 0.016 

m 

Eu 

0.20 

0.0016 

<0.20 

> 0.024 

m 

Dy 

0.20 

0.0007 

<0.20 

> 0.010 

m 

Y 

0.05 

0.0015 

^0.05 

^ 0.090 

m 

Y 

0.10 

0.0015 

<0.10 

> 0.045 

m 

Y 

0.20 

0.0034 

<0.20 

> 0.051 

m 


since our data are determined from fully relaxed configurations, 
which correspond to zero Kelvin. For the proton quantum fluc¬ 
tuations can also play a role. 

6. Conclusions 

Density functional theory (DFT) is used to compute the de¬ 
fect induced strain tensor A for two point defects in barium zir- 
conate: a -\-2 charged oxygen vacancy and a proton interstitial 
(forming a hydroxide ion). The defects have been considered 
both independently and in association with dopant ions. The 
tensor A provides a natural generalisation of the chemical ex¬ 
pansion coefficient to anisotropic deformations. 

We find that the lattice contracts both when the vacancy and 
the proton are inserted, indicating that both the vacancy and the 
hydroxide ion are smaller compared with an oxygen ion. We 
also find that the effect from dopant association is quite small. 
The contraction for the vacancy is considerably larger than for 
the proton and the net effect in hydration, when a vacancy is 
filled and two protons are added, is an expansion, consistent 
with the experimental findings. 
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